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, We discuss the classical motion of a spring of arbitrary mass coupled to two arbitrary massive blocks attached 

O ' at its ends. A general approach to the problem is presented and some general results are obtained. Examples 

* ^ ' for which a simple elastic function can be inferred are discussed and the normal modes and normal frequencies 

obtained. An approximation procedure to the evaluation of the normel frequencies in the case of uniform elastic 

' function and mass density is also discussed. 
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1 Introduction 



The motion of one or two masses linked by a massless spring constrained to move on a straight Une and without 
friction is analysed in several introductory and undergraduate mechanics textbooks, see for example [D, lEI- In 
the case of two arbitrary masses, the two-body problem is solved by a reduction to the problem of a single body 
oscillating with an angular frequency equal to \J Kf- j [i^ where is the spring constant and /i is the reduced mass 
of the system, and the motion of the centre of mass of the system the velocity of which is constant if no external 
forces are present. Moreover, since only the masses make contributions to the kinetic energy and to the total linear 
momentum, the mechanical energy conservation theorem and the linear momentum conservation theorem can be 
applied without much ado. The forces acting on the masses are due to the spring deformation at the extremities 
where the masses are attached to. This is the reason why Newton's third law of motion cannot be directly applied 
to them. We are forced to consider in more detail the mechanism of interaction between the two masses and in 
particular their interaction with the extremities of the spring to which each one of them is attached to. However, 
due to the fact that the spring is massless we can state also that at any given moment of time the sum of those 
forces is zero. Then, in an equivalent way, we can think that the masses move under the action of the force that 
one mass exerts on the other, thereby complying with the third law in such a way that we can ignore the existence 
of the spring. 

The correction to the frequency for the case where one of end of the spring is kept fixed and the mass m of the 
spring, though not zero, is much less than the mass il/ attached to the oscillating end is well known. In this case in 
order to get the angular frequency up to first order we can consider the spring massless and replace the mass of the 
oscillating body by an effective mass that is equal to A/ + m/3, see for example [D, see also |3l and references 
therein. 

In this paper we will consider a more general situation. We will consider the problem of two arbitrary masses, 
say M\ and il/2, attached to a spring of arbitrary mass m. The effects caused by the undulatory propagation of the 
massive spring deformation along the spring length will be taken into account. Solutions to particular situations 
such as the ones described above will be considered as appropriate limits of a less particular solution. We believe 
that the approach we take here may be of some pedagogical value for advanced students and instructors as well. 



2 The equations of motion for the massive spring and attached masses 

We begin by estabUshing the equation of motion for the massive spring along a single spatial dimension. In order 
to do so we introduce an auxiliary parameter x that will help us to describe the properties of the spring such as 
for example its tension or its density at a given point. With this end in mind let us consider the spring in a non- 
deformed condition and denote by ^ its natural length. Now we define a one-to-one correspondence between the 
spring viewed as a one-dimensional smooth matter distribution and the closed interval [0, t\ in such a way that 
a; = corresponds to the left end of the spring and x ^ l\a its right end, see Figure ^ To an arbitrary point P on 
the spring there corresponds a point x E [0, £]. The parameter x must not be viewed as a regular spatial coordinate. 
This parameter can be thought of if we wish as an internal degree of freedom of the spring and it is not subject 
to the transformations associated with the one-dimensional Galileo group, for instance, non-relativistic boosts or 
translations. Were the string made up of N discrete masses labeled by a discrete index j running from 1 to N, this 
index would play a role analogous to x. We assume that the correspondence established here holds for any state 
of motion of the spring, exactly as in the case of the discrete model. Now let it be an inertial reference frame S 
and a suitable coordinate system and let us suppose that the spring moves along the w-axis such that the position 
of a point of the spring with respect to S is given by the function u (x, t). Figure^ The tension T at a point of the 
spring is given by |4l 

T(x,t) = «(.)(^-l) (1) 
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where k (x) is the elastic function of the spring which on physical grounds we suppose to be always positive, that 
is k{x) > for any x e [0, £]. In this way, at a given point the force that the right portion of the spring exerts on 
the left portion wiil heT {x,t) and conversely the force that the left portion of the spring exerts on the right portion 
will be — T {x, t). Consider now an element of the spring determined by a; e a; + dx. The resultant force acting on 
this element is 



dF{x,t) = T(x + dT,t) - T(x,t) 
d_ 
dx 



dx 



dx 



(2) 



If p [x) is the linear mass density of the spring then upon applying Newton's second law of motion to the element 
of mass p [x) dx we obtain 



d_ 

dx 



K (x) 



du {x, t) 
dx 



p{x) 



d'^u [x, t) 
dt^ ■ 



Equation Q controls the motion of the spring. It can be simplified by introducing the variable 

^ {x, t) ~ u (x, t) — X, 
then the equation of motion of the spring becomes 



d_ 

dx 



K (x) 



d^ ix,t) 



dx 



= P [x) 



dt^ ' 



(3) 



(4) 



(5) 



A word of caution: though ^ {x, t) is related to the deformation of spring it does note represent this deformation 
directly. 



O 

h- 



P(0) P{x) P{i) 

P{x) 

> 

u {x, t) 



Figure 1: The motion of a point of the spring with respect to an inertial frame is described by the coordinate 
u (x, t). Given a point P of the spring we associated with it the parameter x. This association is independent of 
the dynamical state of the spring. 



Notice that the equation of motion of the massive spring in the form given by equations Q or (|5} is invariant 
under the galilean transformations. In fact, if we go from the inertial system S to the inertial system S ' that moves 
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with velocity V with respect to S the following evident relations hold 



u' {x,t) = u{x,t) + Vt (6) 
du' {x,t) du{x,t) 



dx dx 
d'^u' {x, t) _ d'^u{x,t) 
dt^ ~ dt^ ■ 



(7) 
(8) 



In this sense the gaUlean invariance of the equation of motion of the spring is manifest in accordance with the fact 
that this equation derives from a straightforward application of the principles of newtonian mechanics. 

Let us now consider the coupled masses. Let us model them by means of two point particles one of mass Mi 
coupled to the left end {x = 0) of the spring and the other of mass M2 coupled to the right end {x = t). Making 
use of equation (|5jl and Newton's second and third laws we can write the equations of motion of the masses as 

The complete solution of these equations and of equation Q demands that we prescribe the initial conditions 

5(a;,t = 0) = ip{x)-x (11) 
d£, (x, t = 0) 



dt 



^{x) (12) 



where Lp {x) := u (x, 0) and ij} (x) := du (a;, t = 0) /dt describe the initial position and velocity of the points of 
the spring. 

Our aim is to obtain a general solution ^ (x, t) - or u (x, t) - to the problem and therefore describe an arbitrary 
state of motion of the system i.e., the two blocks plus the massive string. 



3 General solution of the equations of motion 



We begin by solving equation (|5jl by the method of separation of variables, that is, we look for a solution of the 
form 

ax,t)=X{x)T{t) (13) 

that satisfy also the boundary conditions given by equations and J10> . Taking equation J13l l into (|5Jl and 
introducing the separation constant — A we have 



<fT{t) 



+ XT{t) = 



d_ 

dx 



k{x) 



dX{x) 



dx 



\p{x)X{x) 







(14) 
(15) 



Equation (|5Jl imposes a boundary condition on equation \\5\ to obtain it we take equation \Yi\ into equation (|5Jl 
and write 



MiX (0) 



dt^ 



dx 



(16) 
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and taking into account equation il4\ we obtain 

k(0)X'(0) = -AMiX(O) (17) 
In the same way, taking equation il3\ into dlOt and combining with equation il4\ we obtain the condition 

K {£) X' {£) ^ XM2X {£) (18) 

Let us show now that the eigenvalue A cannot assume negative values. Suppose that we have an eigenfunction 
X{x) corresponding to a particular eigenvalue A. Consider the following identity which can be derived after an 
integration by parts and use of equations ( I17> . ( I18> . and il5\ 



K{x)X'^{x)dx = A 







MiX^ (0) + A'hX^ (£)) + / p{x) X^ (x) dx 



(19) 



Since the left-hand side is always non-negative and the bracket on the right-hand side is always positive we con- 
clude that A is non-negative. 

The null eigenvalue is physically acceptable and has a special meaning. The reason is that A = is common 
to all springs regardless of their mass density, elastic function and the masses of the attached particles at their 
extremities. Notice that the corresponding eigenfunction (the zero mode) can be obtained in a general form. This 
particular mode is not associated with an oscilatory motion of the spring. In fact, for A = the temporal function 
has the form 

r(t) = aoi + /3o, (20) 
where ao and Pq are constants. On the other hand, equation J15l l for A = yields 

dX{x) 7 

where 7 is an integration constant. Boundary conditions as expressed by equations (|9jl and ilO\ . or equivalently 
equations ( I17> and dlSl l. demand 7 = 0, so that X {x) = b = constant. The eigenfunction corresponding to this 
eigenvalue is then 

^oix,t) = b{aot + /3o) 

= xa + Vt (22) 

where we have introduced the new constants xq and V. It is clear that this solution corresponds to a uniform 
motion of the entire system (masses plus spring) with a common velocity V. The zero mode motion is related to 
galilean boosts and may be added to any other solution of the problem if questions about galilean invariance are an 
issue. 

Finally, let us consider the case of positive A. Setting X ~ oj'^ for convenience we write the solutions to equation 
(EJas 

r„ (t) = A„ cos {iOr^t + 0„) (23) 

where n is a positive integer, uJn is the ?i-th frequency eigenvalue indexed in crescent order {uJi < UI2 < W3 . . . ) 
and An and 0„ are constants. The n-th eigensolution to equation|5} corresponding to the n-th eigenfrequency is 

U [x, t) = Xn (x) T„ (t) (24) 

These modes represent the oscillatory modes of the system. The general solution can be written as 

00 

e {x, t) = x-o + + ^ X„ {x) r„ (t) (25) 
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Consequently, in terms of the function u{x,t) the general solution will be given by 



u{x,t) = xo + x + Vt + J2 [x] T„ (t) (26) 

n=l 

The next step is the explicit determination of the spectrum of eigenfrequencies uj. This is, however, a hard task 
to perform and in principle it can be accomplished only if we also know explicitly the elastic function k (x). As 
mentioned before, the zero mode is the only mode that does not depend on the form of k [x). 



4 The orthogonality of the eigenfunctions 



Before dealing with concrete examples let us consider a little bit more some of the formal aspects of our problem. 
Equations ( I17> and (I18> can be read as boundary conditions for equation ( I15t therefore only for certain values 
of A there will be solutions to this equation. The reader will recognise immediately that we are dealing with 
a Sturm-Liouville system. Let us consider then two different eigenvalues, say A™ and A„ and their respective 
eigenfunctions (x) and X„ (x). These eigenfunctions satisfy the differential equations 



d_ 

dx 
d_ 
dx 



K {x) 



dXra {x) 



K {x) 



dx 

dXn {x 



dx 



f XmPXrn (x) = 
+ XnpXn (x) = 



(27) 
(28) 



As usual we multiply the first equation by Xn and the second by Xm, subtract one from the other and after simple 
additional manipulations we end up with 



(A„j - Xn)p (x) X„i (x) Xn {x) + — 

dx 



dx 



dx 



= 



(29) 



Integrating this last equation over the domain [0, £] and taking into account the boundary conditions given by 
equations ( I17> and dlSt we obtain after some simplifications 



(Am ^ A„) 



p (x) X^ (x) Xn (x) dx + (0 X„ (/) + MiX,„ (0) X„(0) 



0. 



(30) 



At this point we define a scalar product in the space of functions that will be convenient for our purposes. Let the 
functions / (x) and g (x) be defined in the closed interval [0, 1], then by definition their scalar product is 



(/,<?>= / p{x)f {x)g{x)dx+ M^f (1)9(1) + M,f{Q)g{0). 



(31) 



With this definition for the scalar product we can consider the eigenfunctions corresponding to different eigenvalues 
as an orthonormal set of eigenfunctions, i.e. 



(Xni 7 Xn 



m, n = 0, 1, 2, 



(32) 



By making use of the initial conditions and the above orthonormality condition the determination of the constants 
An and (f>n that appear in equation i23i and therefore in the general solution can be done in a systematic way. For 
the zero mode, for instance, we have 



Xo(x) 



1 



(33) 



6 



5 Conservation laws 



Linear momentum and mechanical energy conservation theorems can be proven under quite general conditions. 
The former depends on the fact that the system is isolated and the latter depends also on the fact that the internal 
forces can be considered as conservatives. Let us consider first the Unear momentum of the system. Our goal will 
be to determine explicitly the contribution of the massive spring to the total linear momentum. 
The linear momentum due to the two blocks is given by 

P„ = „.2!iM + „,My) (34) 

ot ot 

Making use of equations (|9} and (I10> . we can recast the total time derivative of P^i^t^ into the form 

^-<»>^-<')^ 

On the other hand we can integrate equation (|5} over the domain [0, P\ to obtain 



^ / p(a:)^(x,i)d^ = fcW^^-MO)^^ (36) 



Taking this result into equation (I35> we have 



' 'P...+ /V(x)^rf.)=0. (37) 



dt \ 7o dt 

Defining the linear momentum of the spring by 

, du{x,t) f' , , di{x,t) , 
P.,^, = p{x) -^f^ = P{x) dx (38) 

we see that the total linear momentum of the system Pu^t, + P,p,i„, is conserved. The total linear momentum can be 
rewritten in the form 

,^ , r dui (t) , ^ du2(t) , , duix.t) , 

P.,u, = Ml -^jy- + A'h + p{x) dx (39) 

dt dt Jq ot 

where ui (t) = u (0, t) and U2 {t) = u {£, t). In terms of ^ [x, t) we have 

^.,., = M,^il2 + M,^jl2+/'.(x,^^ (40) 

OT dt In ot 



From equation (I25> or i26\ we can rewrite the total linear momentum in the form 



= {Ml + M2 + m) V + y] MiXn (0) + AhXn (i) + p (x) X„ (x) dx f (t) (41) 

n=l L 

This expression can be rewritten in the form 

00 

P.,,. = (Ml + A/2 + m)V + J2 i^m, Xo {x))f{{t) (42) 
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Since Xq and X„ are orthogonal we see that only the zero mode contributes to the total linear momentum 

i^„,„, = (Ah + M2 + Af,p,,„,) V (43) 

From this result we see that the constant V is the velocity of the centre of mass of the system, as expected. 
We now consider the mechanical energy of the system. The kinetic energy of the blocks is given by 



du{Q,t) 



1 



du {£, t) 
dt 



(44) 



which evidently is not per se a conserved quantity because the blocks exchange energy with the spring. It follows 
that in order to have conservation of the mechanical energy it is mandatory that any variation of the kinetic energy 
of the blocks be compensated by a variation of the energy of the spring, kinetic, potential or both. Keeping this in 
mind we derive ll44> with respect to the time to obtain 



- m + T,) _ Mi^^ + Ah^^ ^ 

Combining this result with equations (|9j and (I10> we can eliminate the masses of the blocks and write 



(45) 



k{0) 



du (0, t) 

at 



du (0, t) 
dx 



du {£, t) 
dt 



du {£, t) 
dx 



(46) 



We can recast this equation into a more useful form if we first multiply equation ^ by du {x, t) / dt to obtain 



p{x) 



du {x, t) d^u (cc, t) 



d_ 

dx 



K (x) 



/ du {x, t) 
\ dx 



- 1 



du [x, t) 



dt dt^ 

Then integrating this result over the interval [0, £] we will have 

'du {i,t) 



dt 



1 , . d f du{x,t) 

-r^'^^di 



dx 



- 1 



(47) 



dx 



1 



du{l,t) .^.duiO,t)du{0,t) 
K (0) ■ 



dt 



dx 



dt 



+ 



d 
di 



-K (x) 



Taking this last equation into equation (I46> it follows after one integration more that 



E ■ 



du {x, t) 
di 



dx 



-K (x) 



du (x, t) 
dx 



du {x, t) 
dt 

du (x, t) 
dx 



1 dx 



dx 



(48) 



(49) 



This equation expresses the conservation of the total mechanical energy of the system. The first three terms on the 
rhs of equation J49> represent the kinetic energy of the blocks and of the massive spring, the last term represents 
the potential energy of the spring. We can rewrite equation W)\ in terms of the ^ (x, t) 



E = Ti+T2 + 



ix,t) 
dt 



dx 



-k{x) 



ix,t) 
dx 



dx 



(50) 



which turns out to be more useful in some applications. 
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Proceeding as in the case of the total Hnear momentum we can write the total energy in terms of the general 
solution. The result is 



oo 



MiXn (0) + M2X„ {£)+ p {x) Xn [X) dx 



E 

n,m— 1 
/•oo 



MlX„ (0) X,n (0) + AhXr, {£) X,n {£) + R {x) X„ {x) X,„ [x) dx 



Tn (t) 
Tn (t) tn {t) 



/ k{x) ^ X:,{x)X:,,{x)T„ {t)T,n{t) dx 
Jo , 



(51) 



The last term representing the potential energy of the spring can be suitable rewritten with help of the following 
identity 



K {x) Xj[ (x) X^ {x) dx — X„ 



MiX„ (0) X^ (0) + A'hXn {£) X^ {£)+ p [x) X„ [x) X^ {x) dx 



(52) 



which can be easily proven. The final result is 



MiX„ (0) + M2X„ {£)+ p (x) Xn {x) dx 



E 

n,m— 1 



MiX„ (0) X^ (0) + AhXn {£) X,n {£)+ p {x) X„ (x) X„, [x) dx 



Tn{t) 
Tn (t) tn {t) 



poo 

/ n{x) (Xn (x) Xm {x)) T„ (t) T„, (t) dx 
Jo „ 



(53) 



^ n.7n—l 

Taking into account the orthonormality relation we can recast the total energy in the more illuminating form 

£^,„ua = 2 (A/i + m2 + M,,„„J V^ + Y, [tl (0 + t^] (54) 



n=l 



or 



^ 00 
E,„,^ = - (Ml +1112 + A/,p,„,) V^ + Y («« + 



(55) 



This last result shows that the total energy of the system can be decomposed into a sum of energies each one 
associated with a normal mode. Moreover, we can see that in order to excite two or more frequencies of comparable 
amplitudes it is necessary to supply the mode with the highest frequency with a greater amount of external energy. 



6 Solution for (x) = 

When the spring is massless the motion of the two blocks is easily obtained by reducing the two-body problem 
to the the motion of a single effective body about a centre of force 13 . Here we try to obtain those solutions by 
making use of equations ( I14> and ( I15> . 
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Firstly, notice that taking p (x) = 0, does not eliminate the possibility of having eigenvalues different from 
zero. It only means that the spatial eigenvalue equation is 



d 

dx 



dX{x) 



the solution of which is 



X{x)=C 

The position of a point of the spring is then given by 



u (x, t) = Xq + X + Vt + 

Defining the usual spring constant by 



C 



dx 
dx' 

K {x') 

' dx' 
dx 





= 



XiO) 



XiO) 



A cos {ujt + (j)) 



with the help of equation \51\ we obtain 

C = K, [X {£) - X (0)] 

Making use of the boundary conditions in the form given by equations ( I17> and ( I18> we have 

KeX (£) - (/Ve - w^A/i) X{0) = 
{Ke - cj^A/a) X (i) - KeX (0) = 



(56) 



(57) 



(58) 



(59) 



(60) 



(61) 
(62) 



In order to have a non-trivial solution the determinant associated with this linear system must be zero, that is 

M1M2 uj^ - Ke (Ml + M2) uj"^ = (63) 
It follows that the allowed eigenfrequency is as expected given by 




where 



M1M2 



(64) 



(65) 



M2 + Ml 

is the reduced mass of the system block 1 plus block 2. This eigenfrequency and the zero mode frequency tJo = 
are the only allowed frequencies of the system when the spring is massless. Making use of equation of our 
definition of the scalar product we can calculate easily the constant C that normahses the eigenfunction. The result 
is 

C=^ (66) 



10 



7 Solution for k{x) and p{x) uniform 



We now turn our attention to an important special case. When the elastic function of the spring and its density are 
uniform it is possible to solve analytically the equation of motion, i.e. the wave equation that describe the system 
and interpret clearly the solutions. Defining 

k/p (67) 
and with A := cj^ to be in accordance with the standard notation, equations ilAl and ( I15l l read 

^+^'T = (68) 

^ + 9'^ = (69) 

where 

LO^qv. (70) 
The boundary conditions, equations and or ( I17> and ( I18> applied to this particular case lead to 

MiX (0) q^v^ + kX' (0) = (71) 
M2X (i) q^v^ - kX' {£) = (72) 

The general solution for the spatial part is 

X (x) ^ Acqs qx + B sin qx (73) 

The allowed eigenvalues are determined by the linear algebraic system 

Miqv^A + kB = (74) 
[M2qv'^ cosq£ + Ksmqe) A+ (M2qv'^ sinqi - Kcosqi) B = (75) 

whose characteristic equation is 

tan q£ = — — ^ (76) 

2 _ P 

Ml + M2 

To illustrate the discussion let us consider the situation for which the density of the spring is very small. In this 
case the mass of the spring can be neglected. Making the necessary approximations to equation (17 6> we obtain 

'^-/^ '''' 
The angular frequency is given by equation (I70> and in this case it leads to 

(78) 

As expected the last result has a non-trivial limit as the mass density tends to zero. Notice that k/£ can be identified 
with the usual elastic constant of the spring. Notice also that the speed of the wave does depend on the density of 
the spring and tends to infinity as the mass density tends to zero. It is precisely this fact that in this approximation 
makes possible to replace the real forces by forces between the two point masses obeying Newton's action and 
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reaction principle discussed in the introduction. To investigate the next order correction to the angular frequency 
we consider for simplicity the case where one of the point masses, say Mi, is infinite and the total mass of the 
point particles Mi + M2 is also infinite. This situation corresponds to the case where one of the extremities of the 
spring is fixed to a wall. Adding one more term to the expansion of tan q£ in ( I76l l we obtain the following quartic 
equation for q 

]<l' + - - (79) 



3^ i^^ M2£3 



whose physical solution is given by 



Consequently, we will have 



A little bit more of simple algebra allows to write 



a well known result, see for example (T]. |3). 




i (M2 + im) 



(80) 



(81) 



(82) 



8 The angular eigenfrequencies 

In order to investigate a general solution of equation ( I76> let us define the variable z = q£ and write the character- 
istic equation ( I76> in the form 

fiz m 
m Mz 



cot z 



(83) 



where M := Mi + M2- In Figure|2we plot the Ihs and the rhs of equation j83t separately for representative 
values of /i, M and m. The solutions of the charcteristic equation are determined by the intersection points. It is 
easily seen that there is an infinite number of eigenfrequencies, one in each open interval (nvr, (n + l)7r), where 
n is a non-negative integer. The lowest eigenfrequency lies in the interval (0, tt). The lowest eigenfrequency is 
the only one that remains finite when the mass of the spring tends to zero. All other eigenfrequencies tend to 
infinite and this means that they are increasingly harder to excite. For n ^ the highest eigenfrequencies can be 
approximately described by the simple formula Zn — n-K- Then we can write 



(84) 



In order to obtain an analytical approximate solution for the eigenfrequencies we solve equation ( I83t for m to 
obtain 

m = ^ — ( -M ± a/ (a/2 + 4 (tan^ z] iiM) ) z. (85) 




2 tan z 



y^(A/2 + 4(tan2 z) fiM)^ 



where the plus sign must be used if x G {mr, mr + tt/2) and the minus sign when if 2; G {mr + tt/2), (n + 1)tt. 
Now we define w = y/m and make use of the Biirmann- Lagrange theorem |5| to express the inverse function in 
the series form. The result is 



1 



1 

6^ 



1 

2M 



w 
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1 



360//2 12/zM 8M2 



w 



(86) 



12 



10 
8 
5 
3 

N 
>^ -3 
-5 
-8 
-10 







V 


V 











6 9 
z 



12 



Figure 2: The dashed curve is the graphical representation of the rhs of equation J83t for given m, M\ and M^. 



Consider only the first term of this series. Then it is easily seen that 

1 Pm 



In this case the angular frequency is 
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(87) 



(88) 



Let us consider the first correction to this result which means to take into account the term in w'^ in the inverted 
series. Then it follows that 
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The angular frequency is then 
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We can also express the other eigenfrequencies (n > 0) in a series form by using again the Biirmann-Lagrange 
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Figure 3: The curves show the behaviour of the frequency of the first three modes of the system as a function of 
the mass of the spring. The sum of the attached masses is same for the three curves and in this case one of the 
attached masses is a fixed wall. 



theorem. The result up to the fourth power in the mass of the spring is 




1 



(91) 



Figures |3]|3 and |5] show the behaviour of the first three lowest eigenfrequencies as a function of the mass of the 
spring m for a particular choice of the sum of the attached masses Mi + M2- In Figures |3]and|3one of the attached 
masses is infinite and the other one is finite. This means that one of the ends of the spring is attached to a fixed 
wall. In Figure|5]both attached masses are finite. As physically expected, when the mass of the spring goes to zero 
the higher modes become harder and harder to excite and the fundamental mode tends to a fixed value. 



9 Final remarks 

In this paper we discussed the classical mechanics of a spring of arbitrary mass coupled to two arbitrary massive 
blocks attached at its ends. A general approach to the problem was attempted and some general results such as 
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Figure 4: The curves show the behaviour of the frequency of the first three modes of the system as a function of 
the mass of the spring. One of the attached masses is a fixed wall. 



the conservation of linear momentum and energy were obtained. We have shown also that the physical problem 
leads to an example of a Sturm-Liouville system. The detailed study of this problem is heavily dependent on 
the explicit knowledge of the elastic function k (x). The special case for which the elastic function and the mass 
density are uniform was discussed and an approximation procedure to the evaluation of the normal frequencies was 
put forward and tested. In the limiting case of a massless spring, we have focused our attention on the motion of 
the attached masses Mi and M2, and considered the spring as a way of transmitting the interection between them. 
With the respect to the wave motion of the spirng, we observe that the result \fK^Jp is the velocity of the wave 
only if the velocity of the centre of mass of the system is zero. If this velocity is V with respect to some suitable 
reference frame then accordding to the galilean rule the velocity of pulse propagation will be ^ = ^ Kf.j p. 

At the moment, the study of a possible equivalence between motion in a single mode of the massive spring and 
simple harmonic motion and possible quantisation of the system is under way. 
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Figure 5: The curves show the behaviour of the frequency of the first three modes of the system as a function of 
the mass of the spring. The sum of the attached masses is same for the three curves and in this case both attached 
masses are finite. 
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